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Abstract 



The field-theoretical renormalization group approach is used to estimate the 
universal critical value <?g of renormalized sextic coupling constant for the 
two-dimensional Ising model. A four-loop perturbative expansion for g§ is 
calculated and resummed by means of the Pade-Borel-Leroy technique. Under 
the optimal value of the shift parameter 6, providing the fastest convergence of 
the iteration procedure the estimates g§ = 1.10, g^/g^ 2 = 2.94 are obtained, 
which agree quite well with those deduced recently by Zinn, Lai, and Fisher 
[Phys. Rev. E 54, 1176 (1996)] from the high-temperature expansions. 
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I. INTRODUCTION 



The two-dimensional (2D) Ising model in a zero magnetic field was exactly solved by 
L. Onsager more than 50 years ago. Since then, many attempts were made to develope a 
theory describing the behavior of this model under an external field but, nevertheless, an 
exact analytical description of its magnetized state has not been given up to now. On the 
other hand, for decades the critical behavior of various systems in ordering fields has at- 
tracted constant attention, being of prime interest both for theorists and experimentalists. 
Recently, the free energy (effective action) and, in particular, higher-order renormalized 
coupling constants #2fc for the basic models of phase transitions, have become the target of 
intensive theoretical studies |I|-I4[. These constants are related to the non-linear suscepti- 



bilities X2k, which enter the scaling equation of state and thus play very important role at 
criticality. Along with critical exponents and critical amplitude ratios, they are universal, 
i.e., possess, under T — > T c , numerical values that are not sensitive to the physical nature 
of phase transition, depending only on the system dimensionality and the symmetry of the 
order parameter. 

Calculation of the universal critical values of ge, g%, etc. for three-dimensional Ising 
and 0(n)-symmetric models by a number of analytical and numerical methods showed that 
the field-theoretical renormalization group (RG) approach in fixed dimensions yields most 
accurate numerical estimates for these quantities. It is a consequence of a rapid convergence 
of the iteration schemes originating from renormalized perturbation theory. Indeed, the 
resummation of four- and five-loop RG expansions by means of the Borel-transformation- 
based procedures gave the values for g>|, which differ from each other by less than 0.5% |||| 
while the use of resummed three-loop RG expansion enabled one to achieve an apparent 
accuracy no worse than 1.6% |]7|JT3|. 

The field-theoretical RG approach, being very effective in 3D (see, e.g., Refs. [|T^-[T^,|TB|) is 
also known to be powerful in two dimensions. Properly resummed four-loop RG expansions 
lead to fair numerical estimates for critical exponents of a 2D Ising model |15| and give 



reasonable results for its random counterpart [18,19]. It is natural, therefore, to use the 
field theory for the calculation of renormalized higher-order coupling constants of a 2D Ising 
model. In this paper, the 2D RG expansion for renormalized coupling constant will be 
calculated, and the numerical estimate for its universal critical value g$ will be obtained. 



II. RG EXPANSION FOR SEXTIC COUPLING CONSTANT 



Within the field-theoretical language, the two-dimensional Ising model in the critical 
region is described by Euclidean scalar field theory with the Hamiltonian 



H 



d 2 x 



l -ml^ + i(V<^) 2 + A<^ 4 



(1) 



where a bare mass squared tuq is proportional to T — being the phase transition 

temperature in the absence of the order-parameter fluctuations. Taking fluctuations into 
account results in renormalizations of the mass m — > m, the field = ipRsfZ, and the 
coupling constant A = m 2 Z 4 Z~ 2 g4. Moreover, thermal fluctuations give rise to many-point 



2 



correlations < <p(xi)(p(x2)-..tp(x2k) > and, correspondingly, to higher-order terms in the 
expansion of the free energy in powers of magnetization M: 



F(M) - F(0) = m 



1M 2 (M 2 \ 2 » /M 2N 



(2) 



In the critical region, where fluctuations are so strong that they completely screen out the 
initial (bare) interaction, the behavior of the system becomes universal and dimensionless 
effective couplings g2k approach their asymptotic limits #| fc . 

In order to estimate g$, we will calculate RG expansion for g e and then apply Pade-Borel- 
Leroy resummation technique to get proper numerical results. Accurate enough numerical 
estimates, as is well known, may be extracted only from a sufficiently long RG series. Below 
we will find the expression for g 6 in the four-loop approximation that will be shown to 
provide an interesting amount of fair numerical estimates. 

The method of calculating the RG series that we use here is straightforward. Since 
in two dimensions higher-order bare couplings are believed to be irrelevant in RG sense, 
renormalized perturbative series for g§ can be obtained from a conventional Feynman graph 
expansion of this quantity in terms of the only bare coupling constant - A. In its turn, A 
may be expressed perturbatively as a function of renormalized dimensionless quartic coupling 
constant g 4 . Substituting corresponding power series for A into original expansion, we can 
obtain the RG series for g 6 . As was recently shown |§,0,|J the one-, two-, three- and four- 
loop contributions are formed by 1, 3, 16, and 94 one-particle irreducible Feynman graphs, 
respectively. Their calculation gives: 



36 [XZ 2 
9e = — \ — 2" 



2' 



XZ 2 XZ 2 \ I XZ 

1 - 11.817855 — — + 110.37270 — - - 985.575 — — 
m 2 \ m 1 / \ m l 



(3) 



The perturbative expansion for A emerges directly from the normalizing condition A 
vr^Z^Z^g^ and the well-known expansion for Z 4 



Z A = 1 + - g A + 5.1288114 g\ + 10.511670 g\ + 0{gi). (4) 

7T 



Combining these expressions we obtain 



36 

— gl(l- 3.2234882 g 4 + 14.957539 g\ - 85.7810 gf). (5) 



7T 

This series will be used for estimation of the universal number g§ 



III. RESUMMATION AND NUMERICAL ESTIMATES 

Being a field-theoretical perturbative expansion the series (^j has factorially growing 
coefficients, i.e., it is divergent (asymptotic). Hence, direct substitution of the fixed point 
value g\ into ([5]) would not lead to satisfactory results. To get reasonable numerical estimate 
for gfg, some procedure making this expansion convergent should be applied. As is well 
known, the Borel-Leroy transformation 
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/(*) = £ = I ^e-*F(xt)dt, F(y) = £ tt^W- (6) 

i=0 J Q i=Q V 1 + °/' 

diminishing the coefficients by the factor (i + b)\, can play a role of such a procedure. Since 
the RG series considered turns out to be alternating, the analytical continuation of the 
Borel transform may be then performed by using Pade approximants. With the four-loop 
expansion (|5|) in hand, we can construct three different Pade approximants: [2/1], [1/2], and 
[0/3]. To obtain proper approximation schemes, however, only diagonal [L/L] and near- 
diagonal Pade approximants should be employed p0 |. That's why further we limit ourselves 



to approximants [2/1] and [1/2]. Moreover, the diagonal Pade approximant [1/1] will be also 
dealt with, although this corresponds, in fact, to the usage of the lower-order, three- loop 
RG approximation. 

The algorithm of estimating g$ that we use here is as follows. Since the Taylor expan- 
sion for the free energy contains as coefficients the ratios i?2fe = g2k/g^~ 1 , rather than the 
renormalized coupling constants themselves, 

F(z) - F(0) = ^(j + z 4 + Rez 6 + Rsz s + ...), z 2 = (7) 

we work with the RG series for i? 6 . It is resummed in three different ways based on the 
Borel-Leroy transformation and the Pade approximants just mentioned. The Borel-Leroy 
integral is evaluated as a function of the parameter b under g± = g\. For the fixed point 
coordinate the value, g\ = 0.6125 (Refs. [|6l,pl),p2]]) is adopted, which was extracted from 



lengthy high-temperature expansions and is believed to be the most accurate estimate for 
g\ available today. The optimal value of b, providing the fastest convergence of the iteration 
scheme, is then determined. It is deduced from the condition that the Pade approximants 
employed should give, for b = b opt , the values of i?g, which are as close as possible to each 
other. Finally, the average over three estimates for R* & is found and claimed to be a numerical 
value of this universal ratio. 

The results of our calculations are presented in Table 1. As one can see, for b = 1.24 all 
three working approximants lead to practically identical values of R§. Hence, we conclude 
that for 2D Ising model at criticality 

Rl = 2.94, g* % = 1.10. (8) 

How close to their exact counterparts may these numbers be? To clear up this point, let 
us discuss the sensitivity of numerical estimates given by RG expansion fl5|) to the type of 
resummation. The content of Table 1 implies that, among others, the results given by Pade 
approximant [2/1] turn out to be most strongly dependent on the parameter b. This situation 
resembles that of 3D 0(n)-symmetric model where Pade approximants of [L — 1/1] type for 
the /^-function and critical exponents lead to numerical estimates demonstrating appreciable 
variation with b, while for diagonal and near-diagonal approximants the dependence of the 



results on the shift parameter is practically absent [p^|15|,|17[ . In our case, Pade approximants 
[1/1] and [1/2] may be referred to as generating such "stable" approximations for g£. Since 
for b, varying from to 15 (i.e., for any reasonable b), the magnitude of g$ averaged over 
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these two approximations remains within the segment (1.044, 1.142), it is hardly believed 
that the values (H) can differ from the exact ones by more than 5%. 

Another way we propose to judge how accurate our numerical results are is based on 
the comparison of the values of g$ given by four subsequent RG approximations available. 
While within the one-loop order we get g$ = 2.633, which is obviously very bad estimate, 
taking into account of higher-order RG contributions to improves the situation markedly. 
Indeed, two-, three-, and four- loop RG series, when resummed by means of the Pade-Borel 
technique with use of "most stable" approximants [0/1], [1/1], and [1/2] yield for gl the 
values 0.981, 1.129, and 1.051, respectively. Since this set of numbers demonstrates an 
oscillatory convergence, one may expect that the exact value of renormalized sextic coupling 
constant lies between the higher-order three-loop and four-loop estimates. It means that 
the deviation of numbers (|8|) from the exact values would not exceed 5%. 

Thus, we see that the four-loop RG expansion for g 6 and elaborated approximation 
scheme lead to accurate enough numerical data. On the other hand, the above arguments, 
as is always the case when we deal with diverging series, should be interpreted as suggestive, 
i.e., they would help us to fix only an apparent accuracy. It is of prime importance therefore 
to compare our estimates with those obtained by other methods. Recently, S.-Y. Zinn, S.-N. 
Lai, and M. E. Fisher, analyzing the high temperature series for various 2D Ising lattices, 



found that i?g = 2.943 ± 0.007 ||; an almost identical value was obtained later in Ref. [14 



Our result for R* & is seen to be in a brilliant agreement with this number. Of course, the 
practical coincidence of the lattice and four-loop RG estimates is occasional, and cannot 
be considered as a manifestation of extremely high accuracy of the methods discussed. 
The closeness of these estimates to each other, however, unambiguously demonstrates the 
power of both approaches. Moreover, such a closeness shed light on the role of a singular 
contribution to g$, which can not be found perturbatively: this contribution is seen to be 
numerically small. 

It is instructive also to compare our results with those given by another field-theoretical 
approach - the famous e-expansion. Today, for the Ising systems, only three terms in the 



e-expansion for R 6 are known [23 



Rl = 2e ^1 - ^ e + 0.63795 e 2 ^J . (9) 

Let us apply a simple Pade-Borel procedure to this series as a whole and to the series in 
brackets, and then let e = 2. We find R$ = 3.19 and R$ = 3.12, respectively, i.e., the 
numbers differ from our estimate by less that 9%. Keeping in mind the lack of a small 
parameter, these values of i?g may be referred to as consistent. We believe that the proper 
account for higher-order terms in the e-expansion for R Q will make corresponding numerical 
estimates closer to those extracted from 2D RG and high-temperature series. 



IV. CONCLUSION 

To summarize, we have calculated the four-loop RG expansion for the renormalized sextic 
coupling constant g§ of the two-dimensional Ising model. Resummation of this series by the 
Pade-Borel-Leroy method has lead at criticality, under the optimal value of the parameter 
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b, to the results gl = 1.10, g^/gl 2 = 2.94. Having analyzed the sensitivity of the RG 
estimates to the type of resummation procedure, and the character of their dependence on 
the order of the RG approximation, an apparent accuracy of these values has been argued 
to be no worse than 5%. Comparison of the RG estimates with their counterparts given 
by other approaches has shown that they are in very good agreement with those deduced 
recently for 2D Ising lattices from high-temperature expansions, and consistent as well with 
the result given by the Pade-Borel resummed three- loop e-expansion for i?g. 

Notes added in proof, (i) Apart from Refs. 1-14, the higher order coupling constants 
in D dimensions were studied in Refs. |2"3-|2"5f . We are grateful to Dr. S. Boettcher, who 



brought these papers to our attention, (ii) Very good agreement exists between the first 
number Eq. (8) and the estimate R$ = 2.95 ± 0.03 (Ref. 14) obtained recently by matching 
the e-expansion available with the exact results known for D — 1 and D — 0. It may be 
considered as an extra argument in favor of the high effectiveness of the field theory in the 
problem duscussed. 
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TABLES 



TABLE I. The values of R$ obtained by means of the Pade-Borel-Leroy technique for various b 
within three- loop (approximant [1/1]) and four- loop (approximants [1/2] and [2/1]) RG approxima- 
tions. The estimate for b = 1 in the middle line is absent because corresponding Pade approximant 
turnes out to be spoilt by a positive axis pole. 



b 





1 


1.24 


2 


3 


4 


5 


7 


10 


15 


[1/1] 


2.741 


2.908 


2.937 


3.009 


3.077 


3.125 


3.161 


3.212 


3.258 


3.301 


[1/2] 


2.827 




2.936 


2.877 


2.853 


2.838 


2.828 


2.814 


2.800 


2.787 


[2/1] 


3.270 


2.988 


2.936 


2.800 


2.667 


2.568 


2.491 


2.380 


2.273 


2.171 
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